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Abstract 

What is the time-optimal way of realizing quantum operations? Here, we show how important 
instances of this problem can be related to the study of shortest paths on the surface of a sphere 
under a special metric. Specifically, we provide an efficient synthesis of a controlled NOT (CNOT) 
gate between qubits (spins |) coupled indirectly via Ising-type couplings to a third spin. Our im- 
plementation of the CNOT gate is more than twice as fast as conventional approaches. The pulse 
sequences for efficient manipulation of our coupled spin system are obtained by explicit computation 
of geodesies on a sphere under the special metric. These methods are also used for the efficient 
synthesis of indirect couplings and of the Toffoli gate. We provide experimental realizations of the 
presented methods on a linear three-spin chain with Ising couplings. 
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1 Introduction 



Quantum computation promises solution to problems that are hard to solve by classical computers 
12). The efficient construction of quantum circuits that can solve interesting tasks is a fundamental 
challenge in the field. Efficient construction of quantum circuits also reduces decoherence losses in physical 
implementations of quantum algorithms by reducing interaction time with the environment. Therefore, 
finding time-optimal ways to synthesize unitary transformations from available physical resources is a 
problem of both fundamental and practical interest in quantum information processing. It has received 
significant attention, and time-optimal control of two coupled qubits (0J |SJ EJ 0) is now well understood. 
Recently, this problem has also been studied in the context of linearly coupled three-qubit topologies 
©, where significant savings in implementation time of trilinear Hamiltonians were demonstrated over 
conventional methods. However, the complexity of the general problem of time optimal control of multiple 
qubit topologies is only beginning to be appreciated. The scope of these issues extends to broader areas 
of coherent spectroscopy and coherent control, where it is desirable to find time optimal ways to steer 
quantum dynamics between points of interest to minimize losses due to couplings to the environment 

OH. 

One approach of effectively tackling this task is to map the problem of efficient synthesis of unitary 
transformations to geometrical question of finding shortest paths on the group of unitary transformations 
under a modified metric fi | 8 3 ) . The optimal time variation of the Hamiltonian of the quantum system 
that produces the desired transformation is obtained by explicit computation of these geodesies. The 
metric enforces the constraints on the quantum dynamics that arise because only limited Hamiltonians 
can be realized. Such analogies between optimization problems related to steering dynamical systems 
with constraints and geometry have been well explored in areas of control theory 1(111 112(1 and sub- 
Riemannian geometry 1(13(1. In this paper, we study in detail the metric and the geodesies that arise from 
the problem of efficient synthesis of couplings and quantum gates between indirectly coupled qubits in 
quantum information processing. 

Synthesizing interactions between qubits that are indirectly coupled through intermediate qubits is a 
typical scenario in many practical implementations of quantum information processing and coherent 
spectroscopy. Examples include implementing two-qubit gates between distant spins on a spin chain 
l|14l 115(1 or using an electron spin to mediate couplings between two nuclear spins 1(16(1 . Multidimensional 
NMR experiments require synthesis of couplings between indirectly coupled qubits in order to generate 
high resolution spectral information (17). The synthesis of two-qubit gates between indirectly coupled 
qubits tends to be expensive in terms of time for their implementation. This is because such gates are 
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conventionally constructed by concatenating two-qubit operations on directly coupled spins. Lengthy 
implementations lead to relaxation losses and poor fidelity of the gates. In this paper, we develop 
methods for efficiently manipulating indirectly coupled qubits. In particular, we study the problem 
of efficient synthesis of couplings and CNOT gate between indirectly coupled spins and demonstrate 
significant improvement in implementation time over conventional methods. We also show how these 
methods can be used for efficient synthesis of other quantum logic gates like for example, a Toffoli gate 
lfT%l ITT?1 l2"0"| on a three-qubit topology, with qubits 1 and 3 indirectly coupled through coupling to qubit 
2. Experimental implementations of the main ideas are demonstrated for a linear three-spin chain with 
Ising couplings in solution state NMR. 



2 Theory 

A geodesic is the shortest path between two points in a curved space. Under the standard Euclidean 

metric (dx) 2 + (dy) 2 + (dz) 2 , the geodesies connecting any two points in a plane are straight lines and 

geodesies on a sphere are segments of great circles. In this paper, we study the geodesies on a sphere 

under the metric g — ^ dx ^ (note on the unit sphere y 2 = 1 — x 2 — z 2 ). The solid curve in Fig. 1 

depicts the shortest path connecting the north pole (1,0,0) to a point (0, cos <fi, sin </>), under the metric 

g. The dashed curve is the geodesic under the standard metric and represents a segment of a great circle. 

For <fi = ?, the length L of the geodesic under g is 0.62 77T (as opposed to \ under standard metric). We 

call this metric g, the quantum gate design metric. If w represents the complex number w = x + iz, then 

the quantum gate design metric can be written an 

_ \dw\ 2 
9 ~ l-\w\ 2 ' 

I dw 1 2 

It has marked similarity to the Poincare metric n^ruTpp ' m Hyperbolic geometry 1|21[1 . defined on the 
unit disc. 

We show how problems of efficiently steering quantum dynamics of coupled qubits can be mapped to the 
study of shortest paths under the metric g. Computing geodesies under g helps us to develop techniques 
for efficient synthesis of transformations in the 63 dimensional space of (special) unitary operators on 
three coupled qubits. 



We consider a linear Ising chain, consisting of three coupled qubits (spins 1/2) with coupling constants 
J12 = J23 = J and J13 = (see Fig. 2 A). The coupling Hamiltonian between the qubits is given by Q17[) 

H c = 2irJ(I lz I 2z +I 2z I 3z ). 
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Figure 1: The solid curve in Fig. 1 depicts the shortest path connecting the north pole (1,0,0) to a 
point (O,cos0,sin0), under the metric g. The dashed curve is the geodesic under the standard metric 
and represents a segment of a great circle. 

The spin system is controlled by local unitary operations on individual qubits, which we assume take 
negligible time compared to the evolution of couplings (@J). The strength of couplings limits the time 
it takes to synthesize quantum logical gates between coupled qubits. We seek to find the optimal way 
to perform local control on qubits in the presence of evolution of couplings to perform fastest possible 
synthesis of quantum logic gates. For directly coupled qubits, this problem has been solved. For example, 
a CNOT(l,2) gate which inverts spin 2 conditioned on the state of spin 1 requires a minimum of 0.5 J -1 
Here, we focus on the problem of synthesizing the CNOT(l,3) gate between indirectly coupled spins. 
Fig. 2 B shows the energy level diagram for the CNOT(l,3)operation, where the state of qubit 3 is 
inverted if qubit 1 is in a lower energy state, i.e., in state 1. In the literature, various constructions of 
CNOT(l,3) have been considered with durations ranging from 3.5 J -1 to 2.5 J~ 1 (122H . The main result 
of this paper is that the CNOT(l, 3) gate can be realized in only units of time, where L is the length 
of the geodesic under the metric g for (p = ? as depicted in Fig. 1. This is twice as fast as the best known 
conventional approach. The new pulse sequence for CNOT(l,3) is based on the sequence element shown 
in Fig. 3 A. 

The main ideas for discovering efficient new pulse sequence are as follows. The unitary propagator for a 
CNOT gate is 

CNOT(l, 3) = exp{-i| (2I lz I 3x - I u - hx + ^1)}, 

where 1 is the identity operator and Ika is 1/2 times the Pauli-spin operator on qubit k with a £ {x, y, z} 
Q17[) . Since we assume that local operations take negligible time, we consider the synthesis of the unitary 
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Figure 2: Fig. 2 A shows the coupling topology where the first qubit (Ji) and third qubit (7 3 ) are 
coupled only indirectly via the second qubit (I2) with coupling constants J12 = J23 = J- Fig. 2 B shows 
a schematic energy level diagram for the spin system in a static magnetic field in the z direction, which 
determines the quantization axis. 1 and are low and high energy eigenstates of the angular momentum 
operator along the z direction. The effect of the unitary transformation CNOT(l, 3) is shown in the 
figure. 

operator 

U{ z = exp{-i|(/ lz + I 3z + 2I lz I 3z )}, 
which is locally equivalent to CNOT(l,3) but symmetric in qubits 1 and 3. 

For synthesizing Uf 3 , we seek to engineer a time varying Hamiltonian that transforms the various quantum 
states in the same way as Uf 3 does. The unitary transformation U{ z transforms the operators I\ a and 
/3 Q , with a € {x,y} to — 2Ii a I 3z and — 2I lz I 3a respectively Since U( z treats the operators h x ,iy and 
l3x,3 V symmetrically, we seek to construct the propagator Wf 3 by a time varying Hamiltonian that only 
involves the evolution of Hamiltonian H c and single qubit operations on the second spin. The advantage 
of restricting to only these two control actions is that it is then sufficient to engineer a pulse sequence 
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for steering just the initial state l\ x to its target operator —2I\ x I 3z . Other operators in the space 
{Iia, ^3/3, 2/1^73^} are then constrained to evolve to their respective targets (as determined by the action 
of U( 3 ). Our approach can be broken down into the following steps: 

(I) In a first step, the problem of efficient transfer of l\ x to —21\ x I 3z in the 63 dimensional operator 
space of three qubits is reduced to a problem in the six-dimensional operator space S, spanned by the set 
of operators I\ x , 2I\ y I 2z , Zhyhx, ^hyhyhz, ^hyhzhz, and 2I\ x I 3z (The numerical factors of 2 and 
4 simplifies the commutation relations among the operators). The subspace S is the lowest dimensional 
subspace in which the initial state I\ x and the target state —2I lx I 3z are coupled by H c and the single 
qubit operations on the second spin. 

(II) In a second step, the six-dimensional problem is decomposed into two independent (but equivalent) 
four-dimensional time-optimal control problems. 

(III) Finally, it is shown that the solution of these time optimal control problems reduces to computing 
shortest paths on a sphere under the modified metric g. 

In step (I), any operator in the six-dimensional subspace S of the 63-dimensional operator space is 
represented by the coordinates x = (xi,x 2 ,x 3 , x±, x 5 ,x 6 ), where the coordinates are given by the following 
six expectation values: x x = (I lx ), x 2 = (2I ly I 2z ), x 3 = (2I ly I 2x ), x 4 = (4I ly I 2y I 3z ), x b = (<ttiyl 2z hz) , 
and xq = — (2I\ x I 3z ). In the presence of the coupling TL Cl a rotation of the second qubit around the y axes 
(effected by a rf-Hamiltonian TL a — u a (t)TrJI 2y ) couples the first four components x' A = {x\,x 2 ,x%,Xk) 
of the vector x. In the presence of 7i c , a rotation around the x axes (effected by a rf-Hamiltonian 
T~tb = Ub(t)irJI 2x ) mixes the last four components x' B — (xs,X4,X5,xq) of the vector x. Under x or y 
pulses applied to the second qubit in the presence of H c , the equations of motion for xa and xb have the 
same form: 



dx A .B T 
— 7— = 
dt 



( -1 \ 

1 -u a ,b 

U A ,B -1 

\ 1 / 



XA.B- (1) 



Since evolution of xa and xb is equivalent, it motivates the following sequence of transformations that 
treats the two systems symmetrically and steers l\ x (corresponding to xa — (1,0,0,0)) to — 2I lx I 3z 
(corresponding to xb = (0, 0, 0, 1)): 

(a) Transformation from (1, 0, 0, 0) to (0, x 2 , x 3l in subsystem A with y/ x 2 + x' 3 2 = 

(b) Transformation from (0, x' 2 ,x' 3 , -±=) to (0, 0, -j=) in subsystem A (corresponding to (-^=, -j=, 0, 0) 
in subsystem B). 
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(c) Transformation from (-7=, -7= , 0, 0) to x' 3 , x' 2 ,0) in subsystem B. 

(d) Transformation from (-^, x 3 , x' 2 , 0) to (0, 0, 0, 1) in subsystem B. 

The transformations (b) and (c) represent fast y and x rotations of the second spin respectively and 
are realized by hard pulses which take a negligible amount of time: Transformation (b) is accomplished 
by hard 9 y pulse applied to the second qubit, i.e. a pulse with flip angle 6 (where tan# = x' 2 /x' 3 ) and 
phase y. Similarly, transformation (c) can be accomplished by hard 9 X pulse applied to the second qubit. 
Because of the symmetry of the two subsystems A and B, the transformations (a) and (d) are equivalent 
and take the same amount of time with (d) being the time-reversed transformation of (a) and X5 and 
xq replacing x 2 and x\ respectively. Hence the problem of finding the fastest transformations (a)-(d), 
reduces to a time-optimal control problem in the four-dimensional subspace A, asking for the choice of 
UA(t) that achieves transfer (a) in the minimum time. 

In step (III), this optimal control problem is reduced to the shortest path problem on a sphere (under 
metric g) described in the beginning of the section (see Fig. 1). The key ideas are as follows. Let 



x(t) = xi(t), y(t) = ya-IW + x \{t) an d z (t) — Xi(t). Since UA{t) can be made large, we can control the 

X2 (t 

2:3(4 



angle tan 0(f) = arbitrarily fast. With the new variables (x(t),y(t), z(t)), equation reduces to 





X 




-sin 0(f) 
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= 7T J 










dt 
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sin 9{t) 


-cos0(i) 




y 




z 




cos 6(t) 







z 



(2) 



In this system, the goal of achieving the fastest transformation (a) corresponds to finding the optimal 
angle 0(f) such that (1,0,0) is steered to (0, -j=, -^=) in the minimum time. The time of transfer r can 



be written as J Q r \l sin 2 8{t) + cos 2 8{t) dt. Substituting for sin 6{t) and cos 8(t) from J5J), this reduces to 




(i) 2 + (i) 5 



dt, 



where L is the length of the trajectory connecting (1,0,0) to (0, ^7=, ^)- Thus minimizing r amounts 
to computing the geodesic under the metric g. The Euler-Lagrange equations for the geodesic take the 
form = — ( 7r ^) _1 (|) an d ^(^) = (nJ)~ 1 (^). For geodesies originating from (1,0,0), this implies 

^ - (ttJ)- 1 cos(/<); Z - = (ttJ)- 1 sin(/i), 

and 9(t) = ft, for constant / that depends of 4>- Now differentiating the expression —(f) = tan(/£) gives 

£3 



u A {t) 



f 1 zx — xz 
' irj 



The Euler-Lagrange equations imply that along geodesic curves, 



is constant, 



y 1 00 'if 

implying that in Eq. (1), time optimal UA(t) = u is constant. We now simply search numerically for this 
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constant u and the corresponding r that will perform transformation (a) in system Q • From all feasible 
(it, r) pairs, we choose the one with smallest r. This gives r = 0.627J" 1 and = = it = 1.04. 

Evolving for time r with it = 1.04 results in 0(t) = tan -1 = 0.5476. The optimal flip angle 
for the transformations (a) and (d) is therefore 6 — 0.5476. With this, the total unitary operator Uf 3 , 
corresponding to the transformations (a)-(d) can be written in the form 

U{ 3 = U x exp-t-i^z} exp{-i6>/2y}n y (3) 

with Tl x , y — exp{— i7rJT [27 lz J 2z + 21 2z l 3z + ul2x,2y]}- The pulse sequence for the implementation of 
Wf 3 is evident from |JS} and consists of a constant y pulse on spin 2 of amplitude v a = uJ/2 = 0.52J 
for a duration of t = 0.627J -1 followed by a y pulse and then a x pulse each of flip angle 6 — 0.5476 
(corresponding to 31.4°) on spin 2, both of negligible duration. Finally, we apply a constant x pulse on 
spin 2 of duration r = 0.627J -1 and amplitude 0.52 J. The overall duration for the implementation of 
U{ 3 is T = 2r = 1.253J" 1 . 

We now show that U{ 3 is locally equivalent to U{ 3 (and hence to CNOT(l,3)). Therefore CNOT(l,3) can 
also be implemented in a time T = 2t = 1.253 J~ . Let Zi 3, denote the subspace spanned by operators 
{Iiai hp, Iiahfs} with independent a, f3 € {x, y, z} and I 2 denote the space spanned by operators {/2a}- 
It can be explicitly verified that by construction Uf 3 maps X13 to itself and acts identically as U{ 3 . This 
constrains U\ 3 to a local transformation on the space T\$®Ti. We can therefore find local transformations 
U l a oc , Ul oc , and U l c oc such that 

CNOT(l,3) = U l ? c XJ l h oc Ul 3 XJ l ° c . 

These local transformations are readily computed to equal U l ° c = exp{— i^I^y}, Ul° c = exp{i(7r — 
0)I 2y } exp{i(?r - 0)I 2x } exp{if (I lz + I 3z )}, and U l c oc = exp{-i|(7 l2 - I 3z )} exp{-if I 3x }, Furthermore, 
the propagator 

U 13 = exp{-i| (2I lz I 3z )} = U[ oc U 9 13 , (4) 
representing a % rotation under an effective 1-3 Ising coupling can also be generated in time T = 1.253 J" 1 . 

The methods developed above can also be used for efficient construction of trilinear propagators. This 
problem has been studied in detail in (JoJ. The main results of (8} become transparent in terms of 
geodesies under the metric g. Consider the propagator U zvz (k) — exp{ — i2ttk h z l2yl3z}- To synthesize 
this propagator, we again seek to efficiently steer the various states between points as they would transform 
under U zyz (n). For example I\ x is transferred to I\ x cosk, + sinK ^I\ y l2yHz- Consider again the four- 
dimensional operator space A defined exactly as before. The goal now is to steer the initial vector xa 
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from initial state (1,0,0,0) to (cos k, 0, 0, sin«) in Eq. (1). This in the transformed system Eq. (2), 
reduces to steering (x,y, z) from (1,0,0) to (cos /t,0, shift). The geodesies for this problem have already 
been characterized. The length of the geodesic that maps (1,0,0) to (cos «,0, sulk) under the metric g 
takes the form L(n) — n^—^—^- and leads to a minimum time of t(k) of {^J 



r{«) = (5) 

The time optimal UA{t) in Eq. 1 is again constant and the resulting unitary propagator is locally 
equivalent to U ZI)Z (k). Since all trilinear propagators U a p^ are locally equivalent, we get efficient ways of 
synthesizing all these propagators. 



3 Discussion 

We now compare the efficient synthesis of CNOT(l, 3) presented above with some conventional methods. 
The most simple approach (CI) to implement an indirect CNOT gate between spin 1 and 3 involves 
swapping the state of spins 1 and 2 followed by a CNOT(2, 3) and a final swap between 1 and 2, resulting 
in CNOT{l,3) = SWAP(1,2) CNOT{2,3) SWAP(1,2). The minimum time for implementing CNOT 
gate between directly coupled spins 2 and 3 is 0.5J -1 (QJ and each SWAP operation takes 1.5J -1 units 
of time, resulting in a total implementation time of 3.5J -1 (I22H . 

Since we assume that local operations take negligible time, the central step in synthesizing CNOT(l, 3) 
is the construction of the unitary operator U13 = exp{— in Iizkz}- The standard method of synthesizing 
such an operator creates a trilinear propagator U zzy — exp{— i? ^hzhzhy} and uses the following 
identity U13 = exp{— \Hi\U zzy expjiiJi}, and U zzy = exp{— LHjj} exp{— U\ z liy) expji^}, where 
Hi = \2l2 Z I-ix and H2 = \2l2 X Izy Overall, this takes 2.5 J -1 units of time <|22|l . resulting in a realization 
(C2) of the CNOT(l, 3) gate that takes only 71.4% of the time required for the implementation (CI). 
The time to synthesize CNOT(l, 3) can thus be shortened by reducing the time to synthesize the propa- 
gator U zzy . We can use a more efficient synthesis of these trilinear propagators discussed in © to further 
improve the efficiency of synthesis of CNOT(l, 3). Note 

U zzy = exp{i- / 22 }exp{-iiJ 3:r }exp{-i2i73y}exp{Lff3 ;r }, 

where H^ a = j(2Ii z l2 a + 2l2 a Iiy)i with a £ {x,y,z}. This reduces the implementation time oiU zzy to 
J" 1 and of CNOT(l, 3) to 2 J -1 (C3). This time can be even further reduced by using the fact that the 
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Table 1: Duration t c of various implementations of CNOT(l, 3) 



pulse sequence 


Tc/J- 1 


relative duration 


sequence 1 (CI) 


3.5 


100% 


sequence 2 (C2) 


2.5 


71.4% 


sequence 3 (C3) 


2.0 


57.1 % 


sequence 4 (C4) 


1.866 


53.3 % 


sequence 5 (C5) 


1.253 


38.8 % 



shortest time to produce the propagator U zzy is given by v / 3/(2 J) © (see Eq. and uses the identity 

Z4 Z!y = exp{i- I 2z } exp{-iv3- {2h z I 2x + 2-hxhy + —^hz)} exp{-i— I 2z }. (6) 

This implementation then results in a total time of (2 + \/3)/(2J) = 1.866/ J (C4) for the CNOT(l, 3). 
The implementation (C5) proposed here is still significantly shorter than this. The implementation times 
under various strategies are summarized in Tabled 

We now show how efficient implementation of trilinear propagators can also be used for efficient construc- 
tion of other quantum gates like a controlled controlled NOT (Toffoli) gate on spin 3 conditioned on the 
state of spin 1 and 2 for the linear spin chain architecture. The decomposition given in J20J is based on 
four CNOT gates (requiring 0.5 J~ x each) between directly coupled qubits and two CNOT gates between 
indirectly coupled qubits. Hence, using a SWAP-based implementation of the CNOT(l,3) gates (CI), 
each of which requires 3.5 J -1 , the total duration of the Toffoli gate would be 9 J" 1 (Tl). With the 
most efficient implementation of CNOT(l,3) (C5), each of which requires 1.253 J -1 , the decomposition 
l|2"0")l has a total duration of about 4.5 J -1 (T2). The Sleator-Weinfurter construction l(T§|l of the Toffoli 
gate is based on two CNOT operations between directly coupled qubits, two unitary operations which 
are locally equivalent to the evolution of the coupling between directly coupled qubits, each of duration 
0.25 J -1 and one unitary operator which is locally equivalent to = exp(—ij2Ii z l3 Z ). A naive 

approach for synthesizing \/^i3 using SWAP operations has a duration of 3.25 J -1 , resulting in a total 
duration of the Toffili gate of 4.75 J -1 (T3). Based on the optimal synthesis of trilinear propagators 
© VMi3 can be implemented in 4 \i^ = 1-66 J -1 units of time (see Fig. 3 B). The main identity used 
is \ZUi3 = exp(— i\1l2 Z l-i V ) exp(— z j4Ji 2 /2z^3z) exp(i-|2/2z^3y)- This reduces the overall duration of the 
Sleator-Weinfurter construction to 3.16 J -1 (T4). 

Here, we present even shorter implementations of the Toffoli gate, the propagator of which is given by 
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Utoff = exp(-wr(/i Q lf a /J /3 )) where l\ a = (\ - I u ) and I z / = (§ + hz)- Neglecting terms in the 
Hamiltonian that correspond to single spin operations (as these take negligible time to synthesize), the 
effective Hamiltonian for the Toffoli gate is H to ff = \\21\ z l 2z + 2I 2z l3 X + ^Iizkx + ^hzhzhx}- The 
synthesis of ^{2I lz I 2z + 2I 2z I 3x } is obtained by evolution of couplings for (4J) _1 units of time. In @, 
we showed that the optimal synthesis of the trilinear Hamiltonian \^I\ z I 2z l3 X takes units of time 
(also see Eq. |5J). The term exp(— i^2Ii z I^ x ) is locally equivalent to \fU13 = exp(— i j2Ii z l3z) which can 
be synthesized in 4 4^ — 1-66 J^ 1 units of time, as discussed above (see Fig. 3 B). This decomposition 
results in an overall time for a Toffoli gate of ^j^- = 2.573J" 1 (T5). 

Further time savings in the synthesis of the Toffoli gate are achieved by the following construction (see Fig. 

3 C), which also uses the optimal creation of trilinear propagators. Let U\ = exp(—i^(2Ii x I 2x + 2I 2x l3 Z )), 
U 2 = exp{-iZ(4I ly I 2x I 3z )), U 3 = expHf (2I lz I 2y )) and U 4 = expHf {2I u I 2z + 2I 2z I 3z j). Then it 
can be verified that ex.p(—iH to ff) = U\U 2 U 3 UIU4. Note U 2 and U3 commute and the optimal synthesis 
of U 2 as mentioned before takes units of time, while Ui U 3 and U4 take 0.5J -1 , 0.25J -1 and 
0.25J -1 units of time each. The total time for the synthesis is therefore J~ 1 = 2.16J -1 (T6). This 
is more than four times faster than (Tl) and a factor of 1.46 faster than the optimal implementation of 
the Sleator-Weinfurter construction (T4). In the following section, an experimental realization of these 
methods on a linear 3 spin chain with Ising couplings is presented. 

4 Experiments 

In Fig. 3, schematic representations of the pulse sequences based on sub-Riemannian geodesies are shown 
for the efficient implementation of U13 and \/U\ 3 simulating coupling evolution by angles ^ (Fig. 3 A) 
and j (Fig. 3 B) between indirectly coupled qubits. As shown above, it is straight-forward to construct 



CNOT(l,3) from U13, which differ only by local rotations. Based on ^/Uis, it is also possible to construct 

an efficient implementation of the Toffoli gate (Fig. 3 C). For simplicity, in Fig. 3 it is assumed that 
qubits Ii, I 2 , and ^3 are on resonance in their respective rotating frames, coupling constants are positive 
(J12 = J23 = J > and J13 = 0) and hard spin-selective pulses (of negligible duration) are available. 
More realistic pulse sequences which compensate off-resonance effects by refocusing pulses and practical 
pulse sequences are given in the supplementary material. 

For an experimental demonstration of the proposed pulse sequence elements, we chose the amino moiety 
of 15 N acetamide (NH 2 COCH 3 ) dissolved in DMSO-d 6 l|2*5)l. NMR experiments were performed at a 
temperature of 295 K using a Bruker 500 MHz Avance spectrometer. Spins 1\ and -Z3 are represented 
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by the amino H nuclear spins and spin I2 corresponds to the 15 N nuclear spin. The scalar couplings of 
interest are J12 = —87.3 Hz « J23 = —88.8 Hz and J13 = 2.9 Hz. The actual pulse sequences implemented 
on the spectrometer and further experimental details are given in the supplementary material. 
The propagators of the constructed pulse sequences were tested numerically and we also performed a 
large number of experimental tests. For example, Fig. 4 shows a series of simulated and experimental 
1 H spectra of the amino moiety of 15 N acetamide. In the simulations, the experimentally determined 
coupling constants and resonance offsets of the spins were taken into account. The various propagators 
were calculated for the actually implemented pulse sequences (given in the supplementary material) 
neglecting relaxation effects. In the simulated spectra, a line broadening of 3.2 Hz was applied in order 
to facilitate the comparison with the experimental spectra. Starting at thermal equilibrium, the state 

Pa — hx 

can be conveniently prepared by saturating spins I2 and -Z3 and applying a 90° pulse to spin I±. The 
resulting spectrum with an absorptive in-phase signal of spin I\ is shown in Fig. 4 A. 
Application of the propagator U13 = exp{— i-| 2Ii z I^ z } to pa results in the state 

Pb = 2Ii y I 3z . 

The corresponding spectrum shows dispersive signal of spin 1\ in antiphase with respect to spin J3, see 
Fig. 4 B. 

The propagator y/Uis = exp{— ij 2Ii z l3 Z } transforms the prepared state pa into 

1 / 

Pc = -j=(hx + 2Ii y h z ), 

resulting in a superposition of absorptive in-phase and dispersive antiphase signals of spin I\ , see Fig. 4 
C. 

The Toffoli gate applied to pa yields 

PD = ~^j2^ lx + 21y x I% x — 4:Ii x l2 z l3x)- 

Only the first two terms in po give rise to detectable signals. The corresponding spectrum is a super- 
position of an absorptive in-phase signal of spin I\ and an absorptive antiphase signal of spin I\ with 
respect to spin I2, resulting in the spectrum shown in Fig. 4 D. 

The effect of the CNOT(l, 3) gate can be conveniently demonstrated by using a two-dimensional exper- 
iment ( 24). Fig. 6 shows the resulting two-dimensional spectrum of the 15 N multiplet (corresponding to 
spin J 2 ) which reflects the expected transformations of the spin states of Ii and ^3 under the CNOT(l, 3) 
operation. 



12 



5 Conclusion 



In this manuscript, we have shown that problems of efficient synthesis of couplings between indirectly 
coupled qubits can be solved by reducing them to problems in geometry. We have constructed efficient 
ways of synthesizing quantum gates on a linear spin chain with Ising couplings including CNOT and 
Toffoli operations. We showed significant savings in time in implementing these quantum gates over 
state-of-the-art methods. We believe, the mathematical methods presented here will have applications 

I dw \ 2 

to broad areas of quantum information technology. The quantum gate design metric fiq^jp defined on a 
open unit disc in a complex plane could play an interesting role in the subject of quantum information. 
The methods presented are expected to have applications to recent proposals of making nuclear spins 
acting as client qubits 111 Gil share information efficiently via distributed hyperfinc coupling to an electron 
spin i acting as the bus-qubit. Efficient synthesis of couplings between indirectly coupled spins will 
also be very useful in multidimensional NMR applications to correlate the frequencies of spins 1 and 3 
coupled indirectly through spin 2 l|17[l . Recent numerical optimization studies 11251 12611 indicate that the 
gap between conventional and time optimal methods for synthesis of typical quantum circuits, (for e.g. 
quantum Fourier transform) on practical architectures, increases rapidly with the number of qubits. This 
motivates further mathematical developments along the lines of the present work in searching for time 
optimal techniques of manipulating coupled spin systems. In practical quantum computing, this might 
prove to be very important as minimizing decoherence losses by efficient gate synthesis improves the 
fidelity of gates. Since fault tolerent quantum computing protocols require gates to have a fidelity above 
a certain threshold, optimal gate synthesis methods like presented here could prove critical in practical 
quantum computing. The problem of constrained optimization that arises in time optimal synthesis of 
unitary transformations in spin networks is also expected to instigate new ideas and method development 
in fields of optimal control and geometry. 
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6 Supplementary material 



The H and 15 N transmitter frequencies were set on resonance for spins I\ and I2, respectively. The 
frequency difference between spins I\ and 73 was Ai/i3 = 310 Hz. Selective rotations of the 15 N nuclear 
spin I2 were implemented using hard pulses. Spin-selective proton pulses were realized by combinations 
of hard pulses and delays in our experiments. For example, a selective 90° rotation of spin I3 with 
phase y, denoted 90° (J 3 ), is realized by the sequence element 90°(7i,/ 3 ) - (Ai/2) - 180° (7 2 ) - (Ai/2) - 
180°(/2) 90°_ x (/i,/ 3 ), where Ai = 1/(4Aj/i 3 ) and 90° tx (7i, J 3 ) correspond to non-selective (hard) proton 
pulses, acting simultaneously on 2j and 7 3 . Figures 6 A - 8 A show broadband versions of the ideal 
sequence shown in Fig. 3, which are robust with respect to frequency offsets of the spins. Positive 
coupling constants Jyi = J23 = J > (with J13 — 0) and hard spin-selective pulses are assumed. 
Figures 6 B - 8 B show the actual pulse sequences used in the experiments with the 15 N acetamide model 
system. In the experimental pulse sequences, selective H pulses were implemented using hard pulses 
and delays, where Ai = l/(4Az/i 3 ). Furthermore, the sequences were adjusted to take into account that 
in the experimental model system the couplings J12 and J23 are negative. Broadband implementations 
of weak irradiation periods l|23|l are enclosed in brackets and the number m of repetitions was two in all 
experiments. 
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Figure 3: Efficient pulse sequences based on sub-Riemannian geodesies for the implementation of U13 = 
exp{— i^ 2/12/32} (A), \/Ui3 — exp{— ij 2/i z /3 Z } (B), simulating coupling evolution by angles ? (A) and 
j (B) between indirectly coupled qubits, and of a Toffoli gate (C). Qubits I\, I2, and I3 are assumed to be 
on resonance in their respective rotating frames. The unitary operator U13, which is locally equivalent to 
the CNOT(l,3) gate, is synthesized by sequence A in a total time T£, = 2r = 1.253 J -1 . The amplitude 
of the weak pulses (represented by gray boxes) with a duration of t = 0.627 J -1 is v a = uJ/2 = 0.52 J. 
The hard-pulse flip angles are 8 = 31.4° and a = 180° — — 148.6°. Sequence B of total duration 
(4 + \/7)/4 J -1 = 1.66 J~ 1 synthesizes the propagator \JU\3. The amplitude of the weak pulse (gray 
box) with a duration of v^/4 J^ 1 = 0.661 J" 1 is u w = 3J/%/7 = 1.134 J. Pulse sequence C realizes the 
Toffoli gate in a total time (6 + \/7)/4 J -1 = 2.16 J -1 . The sequence is based on the sequence for \JTJyj, 
and a weak pulse with the same amplitude and duration as in sequence B. 
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Figure 4: Simulated (left) and experimental (right) 1 H spectra of the amino moiety of 15 N acetamide 
with J12 — —87.3 Hz, J23 = —88.8 Hz and J13 = 2.9 Hz. Starting from thermal equilibrium, in 
all experiments the state pa — I\x wa s prepared by saturating spins I2 and I3 and applying a 90° 
pulse to spin I\. (A) Spectrum corresponding to pa = Iix, (B) spectrum obtained after applying 
the propagator U13 = exp{— i^ 2Ji z /3 Z } to pa, (C) resulting spectrum after applying the propagator 
-\ZZ7i3 = cxp{— ij 2Ii z l3z} to pa, (D) spectrum after applying the Toffoli gate to pa- 
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Figure 5: Experimental two-dimensional 15 N spectrum of 15 N acetamide representing the effect of the 
CNOT(l,3) gate. Using the relation CNOT(l,3)=exp{-if (I u -7 3z )}exp{-if 7 3x }Ux 3 exp{-if J 33/ }, the 
experimental pulse sequence was based on the implementation of the propagator U13 shown in Fig. 3 A. 
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Figure 6: (A) Broadband version of the ideal C/13 sequence shown in Fig. 3 A, which is robust with respect 
to frequency offsets of the spins. Positive coupling constants J12 = J23 = J > (with J13 = 0) and hard 
spin-selective pulses are assumed. The delay A is r/(4m) and the flip angle S is 2 / KV a r / {Am) = 0.5119/m 
(corresponding to 29.33°/to). (B) Experimentally implemented pulse sequence synthesizing U13 for the 
spin system of 15 N acetamide with J( 1 H, 15 N) w -88 Hz, m = 2, flip angles a = 148.6°, = 31.4° 
and 5 = 14.66° and delay A = 890.2 fis. Narrow and wide bars correspond to 90° and 180° pulses, 
respectively, if no other flip angle is indicated, z-rotations are represented by dashed bars. 
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Figure 7: (A) Broadband version of the ideal \/U\z sequence shown in Fig. 3 B, which is robust with 
respect to frequency offsets of the spins. Positive coupling constants J%2 = J23 = J > (with J13 = 0) and 
hard spin-selective pulses are assumed. The delay A is \/l j (16m J) = 0.1654/ (mJ) and the flip angle S is 
37r/(8m) (corresponding to 67.5°/m). (B) Experimentally implemented pulse sequence synthesizing 
for the spin system of 15 N acetamide with J( 1 H, 15 N) w -88 Hz, exp{-i (tt/2) h z hz) for J( 1 H, 15 N) rj 
-88 Hz with m = 2, 5 = 33.75°, A = V7/(l6m\ J( 1 H, 15 N)|) = 939.5 /is, Ai = l/(4A^i 3 ) = 806.5 fis and 
A 2 = 1/(2| J( 1 H, 15 N)|) = 5.68 ms. 
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Figure 8: Broadband version of the ideal TofFoli sequence shown in Fig. 3 B, which is robust with respect 
to frequency offsets of the spins. Positive coupling constants J12 = J23 = J > (with J13 = 0) and hard 
spin-selective pulses are assumed. (B) Experimentally implemented pulse sequence synthesizing a TofFoli 
gate for the spin system of 15 N acetamide. Delays A, Ai, A2 and the small flip angle 5 are defined in 
Fig. 7. 
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